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We consider the screening of an impurity moment in a d-wave superconductor under the influence 
of a Zeeman magnetic field. Using the Numerical Renormalization Group technique, we investigate 
the resulting pseudogap Kondo problem, in particular the field-induced crossover behavior in the 
vicinity of the zero-field boundary quantum phase transition. The impurity spectral function and 
the resulting changes in the local host density of states are calculated, giving specific predictions for 
high-field scanning tunneling microscopy measurements on impurity-doped cuprates. 



I. INTRODUCTION 

Impurity effects provide valuable information about 
bulk properties of correlated electronic systems such as 
cuprate superconductors. Remarkably, non-magnetic im- 
purities like Zn and Li can induce a local momcnt-|an. jtlip 
neighboring lattice sites at intermediate energiestmt3'BEl. 
On the other hand, Zn is known to act as pair breaker 
and strong scatterer, e.g. in low-temperature transport 
measurements, and at present it is not clear whether this 
behavior arises from bare potential scattering only, or if 
other effects like Kondo physics play a role. 

Local impurity moments - arising either from magnetic 
impurities or induced by non-magnetic defects - interact 
with the low-energy bulk excitations. For cuprates at 
temperatures T below the superconducting T c , these are 
the fermionic quasiparticles of the d-w&ve superconductor 
and collective bosonic spin fluctuations. The coupling to 
the spin fluctuations mainly plays a role at intermediate 
energy scales of the order of the spin gap (an estimate 
is given by the position of the so-called resonance mode 
at 40 meV) - their effect on the impurity dynamics has 
been discussed recently in Ref. ^. 

For smaller energies, the dynamics of the impurity mo- 
ments is governed by the interaction with the fermionic 
Bogoliubov quasiparticles. Importantly, their single- 
particle density of states (DOS) vanishes at the Fermi 
energy, therefore the physics, is quite distinct from the 
usual Kondo effect in meljaiaj Extensive studies of the 
pseudogap Kondo modelutlliJiIlLj, where the host DOS 
follows a power-law near the Fermi level, p(u>) ~ |c<j| r 
(r = 1 for d-wave superconductors), have shown the ex- 
istence of a boundary quantum phase transition between 
a free moment and a screened moment as the Kondo cou- 
pling^ increased. Applied to the cuprates, we have ar- 
guedti that this can lead to a strongly doping-dependent 
Kondo temperature, caused by the doping dependence of 
the size of the <i-wave gap. These theoretical findings are 
consistent with NMR experimental measuring the local 
susceptibility near a non-magnetic impurity. This sus- 
ceptibility has been found to show Curie- Weiss-like be- 
havior, with a Weiss temperature being finite at optimal 
doping but apparently vanishing for underdoped samples 



- this implies the existence of unscreened moments in the 
underdoped regime down to temperatures of order 1 K. 

Interesting impurity signatures have also been ob- 
served in recent scanning tunneling microscopy (STM) 
experiments of the high temperature superconductor 
Bi 2 Sr 2 CaCu 2 8+( 5 (BSCCO)l4H .Studies in the vicin- 
ity of Zn impurities on the Cu sitesL3 have found a large 
peak in the differential conductance at small energies of 
1-2 meV near the Zn site, with a characteristic spatial 
variation of the signal. This peak could represent .a-quasi- 
bound state in a potential scattering modellijll3'EZl, but 
such a state appears at low energies only for a range 
of large, potential values depending uppa.jnicroscopic 
detailaiH. Recently, it has been proposedliJ'O that some 
of the properties of the Zn resonance can be explained 
by the Kondo spin dynamics of the magnetic moment 
induced by Zn. The screening of this moment by the 
Bogoliubov quasiparticles provides a natural low-energy 
scale (the Kondo temperature) explaining the energetic 
location of the peak seen in STM - in this case the po- 
tential scattering model might still apply as an effective 
model at low temperatures. However, not all features of 
the Zn resonance can be understood within the existing 
models, and important ingredients in the theoretical de- 
scription, e.g., the actual tunnelmg path of the electrons 
from the tip to the Cu0 2 planeta, are not well known, 
calling for further experimental clarification. 

In this paper we study the screening of an impurity 
moment embedded in a d-wave superconductor in a Zee- 
man magnetic field; we do neglect the changes in the 
quasiparticle spectrum due to the superflow induced by 
vortices which is appropriate for in-plane magnetic fields 
(parallel to the copper oxide layers) and for impurities 
located in between two vortices. In particular, we calcu- 
late the field-dependent impurity spectral density - based 
on these results, we propose that STM tunneling experi- 
ments in high magnetic fields will be able to distinguish 
between the possible sources of the observed impurity 
STM signal, and thus contribute to clarify the nature of 
impurity scattering in cuprate superconductors. 
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II. MODEL 

We begin by describing our model Hamiltonian, H = 
HBCS+Hi mp , for a superconductor with a single impurity. 
The first term describes the host superconductor, which 
we model by a simple BCS Hamiltonian 



a) 



b) 



BCS 



'bulk 



A k T*] tt k . (1) 



Here ^k = (ckj,c_ k |) is a Nambu spinor at momentum 
k = (k x ,k y ) in standard notation. t x,v,z are Pauli ma- 
trices in particle-hole space, and /i is the chemical po- 
tential. For the kinetic energy, £k, we have first (i), 
second (t'), and third it") neighbor hopping, while we 
assume a d-wave form for the BCS gap function Ak = 
(A /2)(cosfc x cos/^y). The spin quantization axis lies 
in the direction of the (in-plane) Zeeman field /ibuik- 

The Green's function of the conduction electrons de- 
scribed by Hbcs is G°(k, w) = [u - (e k - p)r z - h huik - 
AfcT x ] _1 in matrix notation. The influence of a Zeeman 
magnetic field on a d-wave BCS superconductor has been 
studied in Ref. [h]. For low temperatures and small fields, 
^buik <C Ao, the field-induced change in the gap function 
Ak can be neglected. Thus the main effect of the Zee- 
man field is to split the quasiparticle spectrum by ±/itmik, 
which leads to a finite DOS at the Fermi level. 

We will not consider the coupling of the external field 
to the orbital motion of the electrons - the primary ef- 
fect of the field-induced superflow would be an addi- 
tional shift in the quasiparticle spectrum. Neglecting this 
Doppler effect is approximately justified for fields applied 
parallel to the Cu02 planes of the system: Due to the lay- 
ered crystal structure and strong anisotropy, the orbital 
effect of the magnetic field is much weaker forja-plane 
fields, and vortex effects are likely to be smalfcl Fur- 
thermore, we expect our approximation to apply also for 
fields in c-axis direction, i.e., perpendicular to the planes, 
for impurities located in between two vortices where the 
superflow tends to cancel and hence the Doppler shift is 
small. 

The impurity, located at ro = (0,0), is described by 
an effective model consisting of a scattering potential U 
and a Kondo term for the spin- 1/2 impurity spin S: 



'in,.. — 



JrS ■ sr + S z h 



loc 



C/4(r )c Q (r ). (2) 



where sr = N 



E Q i(k-k')R„t 

conduction band spin operator at site R. As discussed 
in Refs. ^ 10 1^ we will consider two cases: (i) a local 
magnetic moment where Af = {ro}, and (ii) a spatially 
extended magnetic moment (modelling Zn) where J\f is 
a set of sitesJn the neighborhood of r ; guided by NMR 
experiments] we take N to be the four neighbors of ro, 
i.e., J R = J for R - r = (±1, 0), (0, ±1) and J R = 
otherwise. Note that this four-site model describes an in- 
duced moment fluctuating as a single entity, in agreement 
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FIG. 1: Schematic zero-temperature phase diagrams for the 
asymmetric pseudogap Kondo model in a Zeeman field, de- 
duced from NRG. a) Local field h\ oc only, i.e., for ftbuik = 0. 
There is a line of continuous phase transitions at finite fields, 
terminating in the multicritical point at J = J c , h\ oc = 0. 
(The phases LM, ASC, and POL are explained in the text.) 
b) Uniform field h = h\ oc — ftbuik- Here, the asymptotic low- 
T, low-frequency behavior at finite fields is always that of 
a conventional Kondo-screened spin in a magnetic field (SCR 
phase), but there is a sharp crossover separating regions which 
are dominated by ASC and POL behavior at finite energies. 



with spin relaxation experimentsO. The Kondo term in 
the four-site case defines a multichannel Kondo problem; 
however, the low-energy physics in the present situation 
is dominated by a single channel given by the d-wave- 
like linear combination of-,thc conduction electrons on 
the four neighboring sitestlliJ. 

The potential scattering term in (||) can be accounted 
for exactly: 

G(r, r', u) = G°(r - r', uj) + UG°(r - r , u) 
xt 2 [1- UG°((0, 0), wJrT^ro - r '> w )- ( 3 ) 

The Kondo problem we consider in the following is then 
defined by the spin-dependent DOS seen by the impurity. 
For a single-site impurity we have 



Peff.aM = -ImG Q , Q ,(r ,r ,sgn(a)w) / it 



(4) 



with p e ff,c((w) ~ |w| at zero field, whereas the extended 
four-site impurity leads to 

p e G, a (u) = -J- Im X] <<2R<y9R'G QCt (R,R',sgn(a)w) 



R,R' 



which behaves as p c s, a (^>) ~ M 3 at ft-buik = 0. Here, 
</j R = +[-]l forR-ro = (±1, 0)[(0, ±1)] describe the d- 
wave channel of the four-site impurity, G is the ^ Green's 
function (|3|) in Nambu notation, and sgn(a) = ±1 for 



III. BOUNDARY CRITICAL POINTS AND 
CROSSOVER BEHAVIOR 

In the following, we describe our numerical results for 
the above Kondo model in an external field, obtained 
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FIG. 2: NRG results for the local impurity magnetization 
mioc = (S z ), calculated for a four-site impurity in a cuprate 
band structure (t = -0.15 eV, t' = -t/4, t" = t/lW = 
0.04 eV, ft = -0.14 eV) with J c « 0.23 eV for U = E3, as a 
function of the applied field, a) Local field hi oc only, for a J > 
J c with T* » 5 K. At the phase transition the magnetization 
follows a power law, |mi oc — wiioc,c| ~ \h\ oc — hioc,c\ , with 
£ ~ 0.1. b) Uniform field h — hi oc — ftbulk, for a J > J c 
with T* « 25 K. Here, the magnetization evolves smoothly 
through the crossover region of Fig. [j]b. The lower figures 
show zooms of the upper ones. 



with Wilspn's Numerical Renormalization Group (NRG) 
techniqueo. The zero-field pseudogap Kondo problem 
in the presence of particle-hole asymmetry is known to 
show a boundary quantum phase transition between an 
asymmetric strong-coupling (ASC) phase with a screened 
impurity moment and a local moment (LM)jAhase with 
Curie behavior of the impurity susceptibilitytJO. 

We start by considering the case of a local field only, 
hioc 7^ but /ibuik = 0, such that the host DOS is unaf- 
fected by the field. From the analysis of the fixed points 
of the NRG calculation we have found that even in a finite 
local field there is still a phase transition upon varying the 
Kondo coupling. The T — phase diagram deduced from 
NRG is schematically shown in Fig. |l|a - there is a line of 
finite-field transitions ending jin-thc zero-field transition 
point known from earlier worlfiLil which can be viewed as 
a multicritical point. The three stable phases are easily 
understood: The ASC phase at J > J c extends to finite 
field (with weak spin polarization for h\ oc ^ 0), whereas 
the LM phase at J < J c with a free impurity spin exists 
at zero field only. Applying a field at J < J c immediately 
leads to a phase with strong polarization of the impurity 
spin, therefore we term this phase POL. 

All transitions found are continuous phase transitions, 
characterized by a small energy scale T* , which vanishes 
at the transition, and defines the crossover energy above 
which quantum-critical behavior is observed. In the ASC 
phase at zero field, TX{J, h = 0) can be identified with the 
Kondo temperaturecil, i.e., with the binding energy of the 
Kondo singlet. The transition from ASC to POL is char- 
acterized by breaking this singlet by the external field - 
consequently the critical field h\ c,c close to the multi- 



critical point is given by T*(J, h = 0) times a prefactor 
depending on the DOS exponent r. Note that this novel 
field-induced Kondo transition is unique to the pseudogap 
model studied here - in the metallic host case (r = 0) the 
application of a field-leads only to a gradual suppression 
of the Kondo effectS. 

From the above discussion we expect that a uni- 
form Zeeman field ft-buik = h\ oc = h removes the critical 
point, because a finite /ibuik leads to a finite DOS at 
the Fermi level. Our numerical calculations show that 
this is indeed the case - for any finite h the asymptotic 
zero-temperature behavior is that of a conventional (i.e. 
metallic-host) Kondo-screened spin in a magnetic field. 
However, the characteristic energy scale below which this 
behavior occurs is tiny, as it is exponentially small in the 
field-induced DOS at the Fermi level. (The same applies 
to effects arising from the small low-energy DOS induced 
by a finite concentration of impurities, i.e., the so-called 
impurity bandO.) Therefore, at practically accessible 
temperatures the finite-field behavior is determined by 
the physics of the POL and ASC phases described above, 
and both regimes are now separated by a sharp crossover, 
Fig. |b. This crossover field turns out to be somewhat 
smaller than the critical local field in Fig. [lja,; their ratio 
depends on the field-induced DOS at the Fermi level (i.e. 
the DOS exponent r). 

The properties of the phases are reflected in the field- 
induced impurity magnetization, Fig. |[ This quantity is 
finite for any non-zero field, but shows a power-law sin- 
gularity at the transition in the local-field case (Fig. ||a) . 
Note that the behavior of m\ oc near the transition is not 
universal but depends on J - the reason is that the finite- 
field transitions actually represent a line of critical fixed 
points in the RG sense. In the uniform-field case, Fig. ||b, 
the singularity in m\ oc is replaced by a relatively sharp 
crossover. 



IV. IMPURITY T MATRIX AND STM 

Turning to the dynamic properties of the Kondo im- 
purity, we now discuss the conduction electron T matrix, 
T(u>), which is the important input quantity for the calcu- 
lation of the local DOS measured by STM. The T matrix 
spectral density, PT a {u) = — \mT a {uS) / -k , in zero external 
field has been studied in Ref. [ll] using NRG. It was found 
that pr(u) ~ \lu\ t in the stable LM and ASC phases. The 
critical Rpin t at J c is characterized by pr{w) ~ |^|~ r (for 
r < l)liMHl. The crossover from quantum-critical |w|~ r 
behavior at high energies to \u>\ r behavior at low energies 
leads to a characteristic peak in the spectral density on 
one (!) side of the Fermi level at an energy T*(J, h = 0), 
which can be viewed as Kondo peak (note that the prob- 
lem at hand is intrinsically particle- hole asymmetric). 

For finite field, the spin degeneracy is lifted: applying 
a small field to the system with J > J c splits the Kondo 
peak, in a manner roughly similar to the field-induced 
peak splitting in the conventional (metallic) Kondo ef- 
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FIG. 3: Typical NRG results for the T matrix spectral den- 
sity, *^2i a PT a (uj), at fixed J > J„_for different uniform fields 
h, calculated for a local impuritytj. The height of the Kondo 
peak reaches its maximum in the crossover region of Fig. [Tjb 
and rapidly decreases for higher fields. Inset: Normalized 
weight of the Kondo peak, integrated from u = — 10T* to 
10T*. 

fcct0. For /ibuik = the low- frequency limit of pt{u) 
still behaves as |ci;| r ; of course, finite /ibuik will lead to 
a finite value of pt{u = 0). The simple peak-splitting 
picture is modified when the field becomes comparable 
to the zero- field peak energy, i.e., close to the transi- 
tion/crossover lines in Fig. [I]. In the case of a purely local 
field, we observe that the spectral peak for one spin direc- 
tion approaches the Fermi level and becomes sharper as 
\hi oc — h\ oc< c\ is decreased; it diverges at the phase tran- 
sition as whereas the spectrum for the opposite 
spin broadens and remains non-critical. For a uniform 
field, this effect is again smeared out: in the crossover 
region of Fig. 0b the spectral function shows a peak with 
a finite, small width at the Fermi level, see Fig. ||. Re- 
markably, the weight of the Kondo peak (inset of Fig. |^) 
stays nearly constant at small fields (ASC regime), but 
rapidly decreases (POL regime) once the crossover region 
is passed. 

Applying a field to a system with J < J c immedi- 
ately suppresses the spectral function for one spin species 
whereas the other one shows a smaller peak which moves 
away from the Fermi level and gradually looses intensity 
- so there are no drastic effects as suggested by the phase 
diagram. 

Using the T matrix from the NRG methorS it is 
straightfop/ard to determine the local conduction elec- 
tron DOSEil as measured by STM. Fig. || shows a set of 
calculated on-site tunneling spectra for a four-site impu- 
rity in a cuprate band structure, which is supposed to de- 
scribe the moment induced by a Zn impurity in BSCCO. 
For a zero-field peak position of —1.4 meV, the crossover 
region of Fig. |l|b is reached at an in-plane field of ap- 
proximately 20 T. For larger fields, the peak moves to 
positive bias and looses intensity similar to the inset of 
Fig. |H The spatial intensity pattern of the signal is nearly 
field-independent; it is alternating from site to site, i.e., 
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FIG. 4: Field evolution of the tunneling DOS for the four- 
site Kondo impurity, measured on the im»«rity site, without 
taking into account apossible filter effecttj. The host para- 
maters are as in Fig. bl the Kondo coupling is chosen to yield 
a zero-field peak at —1.4 meV (T* = 15 K), the potential 
scattering is U — 0. The results have been broadened corre- 
sponding to a measurement temperature of 4.5 K; the dashed 
line shows the bulk DOS. (The external field B is related to 
h by h = gfi B B/2.) 



the nearest neighbor of the impurity has a very small 
intensity, whereas the second neighbor has again higher 
intensityOj'LilEj. Fig. ^ shows results for zero potential 
scattering U, we note that no qualitative modifications 
occur for small and moderate U values up to the hopping 
energy t, see also Ref. [llj 



V. CONCLUSIONS 

Summarizing, we have studied the effect of a Zeeman 
field on the Kondo screening of impurity moments em- 
bedded in a host with a pseudogap density of states. 
The application of a local magnetic field leads to a novel 
field-induced boundary phase transition, associated with 
breaking up the Kondo singlet - this transition has char- 
acteristic signatures in the local magnetization and the 
impurity spectral function. In the context of Zn impuri- 
ties in cuprate d-wave superconductors we predict a char- 
acteristic evolution of the impurity STM spectrum with 
applied field which can be tested in future experiments. 
The spectra in Fig. |] can be contrasted to the results 
for a pure potential scattering modelEa, where the scat- 
tering peak is simply split by the Zeeman energy, but 
is not modified in its total weight - therefore these ex- 
periments will help to clarify the general role played by 
non-magnetic impurities in cuprates. 
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